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Abstract 
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1. Introduction 



The Tsallis entropy [1], [2] is a nonenxtensive entropic form which has attracted con- 
siderable attention over the last quarter century [2]. Although the scope of its applications 
is currently unknown, there is a large body of circumstantial evidence pointing toward its 
applicability in diverse directions not only in Physics, but also in the biological and social 
sciences as well as in economics and the humanities. This is one reason why a considerable 
amount of effort has been dedicated to understanding not only the scope but also the prop- 
erties and the dynamical foundations of the Tsallis entropy. 

It is probably fair to state that, in spite of all such efforts and some progress made, 
there is still a lot left to be understood as far as the dynamical foundations of the Tsallis 
entropy are concerned [2]. One way to advance toward this end has been by using meso- 
scopics following the examples of the Langevin and the Fokker-Planck equations [3], [4]. In 
this formalism, methods and results of the diffusion equation have traditionally played a 
dominant role. Lessons learned are extended in the formulation and properties of non-linear 
diffusions [5], of equations with fractional derivatives [6], [7] and of their solutions involving 
the q-analogues of the ordinary transcendental functions [2], [8] - [15]. A wide range, al- 
though still unknown, of applicability of such solutions is guaranteed by a relatively recent 
generalization of the Central Limit Theorem [16] - [19]. There has been a recent resurgence 
of interest in formulating and solving such equations [20] - [25] motivated by the properties 
of the Tsallis entropy. 

Probably the most important distinction between the Tsallis and the BGS entropies is 
the different way they deal with independent events. The BGS entropy is additive, whereas 
the Tsallis entropy is not, generically at least. In order make the Tsallis entropy manifestly 
additive [26], [27] introduced a generalized addition. Alternative proposals were considered 
in [28]. A generalized product distributive with respect to the generalised sums of [26], [27] 
was introduced in [29], [30]. Metric and measure properties of the resulting field, induced by 
the Tsallis entropy when compared to its Boltzmann/Gibbs/Shannon (BGS) counterparts, 
were subsequently examined in [31] - [34]. Despite the successes of the structures of [26], 
[27], [29], [30], one can be left wondering to what extent such generalisations are universal, 
or what parts of the composition property of the Tsallis entropy they really encode. In 
other words, are there any properties which are weakly dependent, or even independent, on 
the underlying field, if that field's operations encode the composition property of the Tsallis 
entropy? 
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We address this question in the present work by combine the approaches of the two pre- 
vious paragraphs. In Section 2, we use one of the fields introduced in [28] as an alternative 
to that of [26], [27, [29], [30] to construct the diffusion equation on it. In section 3, we 
present the full solution to this diffusion equation and we compare it to the solution of the 
ordinary diffusion equation and to solutions of the porous medium equation. We determine 
and compare their asymptotic behaviours and comment on the common features of such 
solutions. This in an attempt to determine properties that seem to persist regardless of the 
base field on which the respective equations are formulated. Such a property appears to be 
the underlying hyperbolicity which is induced by the composition of the Tsallis entropy as 
seen in Section 4, where some general conclusions and speculations are presented. 

2. A DIFFUSION EQUATION AND ASYMPTOTICS OVER 1Z 

The Tsallis entropy [1], [2] is in reality a single parameter family of functionals S q 
parametrized by gGl. Consider a system having a Riemannian manifold M as its sample 
space, endowed with a probability density function p : M — > [0, 1] which is absolutely 
continuous with respect to the Riemannian volume cLvoIm- Most frequently M stands for 
the configuration or the phase space of the system. The Tsallis entropy of p is given by 



where ks is Boltzmann's constant that will be set to one from now on, for simplicity. 
Consider two subsets A, B C M. They are independent if the corresponding probabilities 
Pa, Pb obey the composition law 



where A* B indicates the system resulting from the interaction of A and B. It is immediate 
to check that for such A and B the Tsallis entropy has the composition law 




(1) 



Pa*b = Pa - Pb 



(2) 



S q [pA* B ] = S q [p A ] + S q [p B ] + (l-q) S q [p A ] S q [p B ] 



(3) 



This lead to the definitions [26], [27] of the generalized sum 




(4) 
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and of the generalized product 

x® q y={x 1 - q + y l - q -l} 1 ^, x>0, y>0 (5) 
Unfortunately these two operations do not satisfy the usual distributivity property, since 

x <g) g (y ® q z) ^ (x ® q y) ® q (x <S> q z) (6) 

The question [2] of defining a generalised product () q which is distributive with respect 
to the generalised sum (B q (4) was answered recently [29], [30] where such a product was 
provided in closed form in [29] by 

log[l + (l-qr)3;]-log[l + (l-9r)»] 

(2 - q ) ig^g _ 1 

xO q y = — q (7) 

This gave rise [30] to the field deformation R 9 of R whose metric properties were used to 
reach [30] - [34] some model independent results pertaining to the Tsallis entropy. 

We observe that the generalised sum (4) is similar to the ordinary sum when \x\ ~ 
1, \y\ ~ 1, but it is more akin to the ordinary multiplication when |x| ^> 1, \y\ ^> 1. 
Motivated by this observation and the requirement for distributivity, we proposed in [28] a 
reversal of roles of the generalised addition and multiplication. We obtained, as a result, two 
fields indicated by TZ\ and 7Z 2 in [28]. In the present paper it is sufficient to just use one 
of them for our arguments, let's say TZi, which we will call 1Z from now on. We will keep 
the notation of [28] in the sequel, so the non-extensive parameter will be indicated by k and 
the BGS entropy will be recovered in the limit k — > 0. Let x,y e 1Z. As a reminder, the 

k 

addition © on 1Z was defined [28] by the same law as in (5), namely 

x © y = (x k + y k - l)i (8) 

It turned out [28], that (8) is associative, commutative, with neutral element 1 and the 

k 

opposite of x denoted by was found to be 

Qx = (2 - x k )^ 

k k k 

The subtraction was, as a result, defined [28] by x y = x © (© y) which resulted in 

x e y = (x k - y k + l)s (9) 
The multiplication of 1Z was defined [28] by 

^JW±V+(W)} ! (10) 
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It turned out [28] to be associative, commutative, with identity element (k+1) & and inverse 
element, for x G R \ {1} and for a generic value of k 



k f A; 2 

x = <U + 



x k — 1 



fc k k 

The division was defined [28] as x y = x (0 y) and turned out to be 



i 



An explicit isomorphism u; : TZ — > R is provided by the k- logarithm In*. (a;) [2], namely 



> • y = i 1 + k — - )■ ill: 



w(x) = ^-Jl (12) 

It may be one some interest to compare (12) to the isomorphism r q between R and R g 
which was given in [30] by 

, , (2 - q) x - 1 
r q (x) = 1 (13) 

We observe that in the present case the map / follows a power-law (polynomial) whereas 
r q is exponential. Hence the corresponding fields 1Z and R g are substantially different 
from each other. This difference can be ascribed to the relative reversal of roles of the gen- 
eralised additions in 1Z and R 9 . In TZ (8) reduces to the addition of the Banach space 
l h , for k > 1, whereas in R g (4) is more akin to multiplication for large values of x and y. 

Among all rational powers of x G TZ we will only need the square root ^/x, x > 1 in 
the sequel. It is defined by demanding 

^0^x = x, l<xe7Z (14) 
By using (8), (10) we can easily see that ^/x is related to the ordinary square root by 

^ = {l + y /k(x k -l)y (15) 
We defined [28] the derivative of / : TZ — > TZ as 

= lim {f(y) f(x)} {y x} 

y^tx 



which turned out to be 



D 9 f(x) = {l + ^T^[/(-)] fc }" (16) 
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k k 

The derivative turned out to be linear with respect to ©, ® and also obey Leibniz's 
rule with respect to these two operations. Let 1Z + = {x e 7Z, x > 1}. The exponential 
function [28] exp^ : 1Z — > 7Z + was defined as an appropriately normalized eigenfunction of 
.0® and turned out to be 



exp @ (x) = {1 + ke— j (17) 

The normalization that we used in (17) to specify the integration constant in the eigenvalue 
equation was exp (g) (l) = (1 + A;)^, in complete analogy with the ordinary exponential 
function e x : R — > The integral in 7Z was operationally defined [28] by 

demanding 

D® / f(x) <g> d^x = f(x) 
J® 

k 

where d^x = lim^-^ y Q x. This integral is given by 

f(x)®d®x = jl + Jlifix^-ljx^dxy (18) 

k k 

and, like the derivative D®, it is linear with respect to ©, ® as it should be. Using (16), 
the partial derivatives d^t of functions on 1Z can be defined by 

d&f(x,t) = {l + j^^lfMfY (19) 



and 



d®J(x,t) = i l + ^T^[fM} k Y (20) 



With these definitions, the simplest parabolic equation of the diffusion type that can be 
written on 1Z x 7Z + , is 

k 

d® t f(x,t) = D <g> d (g)X d (g)X f(x, t) (21) 

where D > 1 stands for the diffusion constant. To determine its fundamental solution we 
need to use an appropriate initial condition. As is done for the ordinary diffusion equation, 
we will use a unit source in the initial condition and then proceed to the general solution 
by linearity. The Dirac delta function in 1Z is defined, in complete analogy to R, by 
demanding 

r k i 
5®(x) = 1, l^xeTZ and / 5® ® d®x = (l + fc)* (22) 

J ® 
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which, by using (18), reduces to 

+00 

{1+ J [(8®(x)) k - ljx^dx} 1 * = (fc + 1)* (23) 

—00 

This can be expressed in terms of the Dirac delta function on R by writing 

+00 +00 

j [{5®{x)t ~ Ax^dx = k J 5{x)dx (24) 

—00 —00 

which gives 

[(5®(x)) k - I]**- 1 = kS(x) (25) 



[8*(x)] = <± + k S M\ (26) 



which eventually results in 

= { 1 + k x 

To find the fundamental solution of (21), we impose the initial condition 

f(x a ,t=l) = 5®(x ) (27) 

We find, in complete analogy with the solution to the ordinary diffusion equation in Mxl, 
due to the isomorphism (12) that the solution f(x, t) of the initial value problem of (21) 
and (27) is 



k ® / k k k k f k k k k k 1 

f(x, t)=(Z)\jA®'K®D®t ® exp^ <^ e [x ® x] [4 <g> D ® t] > (28) 

By analogy with the corresponding properties of the fundamental solution on R, we find 
that 

k k k 

{x®x) k = 2® D ®t (29) 

where 

k r k k k 

(x (g> x)k ■= fx (8) x (8) f(x, t) <g> d^x (30) 
J® 

which, upon substitution of (28) in (30), gives 

k 

(x ®x) k ~t (31) 

for t — > oo. If this result is pulled back to R by the inverse of w given in (12), then we 
see that the solution clearly has a diffusive behaviour something which is not unexpected in 
view of the algebraic identification of 1Z and R provided by w. It is far more interesting 
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to see how the asymptotic behaviour of the solution appears from the viewpoint of ' 
be easily checked that, up to an unimportant, for present purposes, constant 



f(x, t) ~ { 1 + k exp 
which for large x and large t simplifies to 

f(x, t) ~ { 1 + k exp 



x 



Ik 



2x k + 1 



((4 J> D) k - l)t k - (4 4 D) k + 1 



It can 



(32) 



x 



2k 



(33) 



((4 <g> D) k - l)t k _ 

The asymptotic solution to the initial value problem is clearly a not diffusion. One can 
observe that the k-th power of f(x,t) behaves like stretched exponential as a function of 
x and t with an additional constant term providing the asymptotically dominant, even if 
trivial, behaviour. This dominant behavior of f k (x,t), is to zeroth order similar to the 
asymptotic temporal behaviour of the solutions to the heat equation on a compact manifold 
[35]. The characteristic exponent determining the rate of approach to the steady state is k 
and the corresponding "generalized diffusion" coefficient is 

\AD) k 



Z>= ±[(4 <£>£>)*-!]; 



D k - 4 fc + 1 
~k 



(34) 



It is may be of some interest to observe that V explicitly depends on A; a fact that 
becomes more relevant particularly when systems, and the resulting fields, described by 
different values of k, interact with each other. 



3. General solution and implications 

In this section we intend to obtain explicitly a formal general solution of (21). We proceed 
by using separation of variables. It might be also of some interest to use a variation of the 
recently introduced q- Fourier transform [17], [36] - [39] and its inverse to that end. We speak 
about a formal solution of (21) since we will not specify any particular initial conditions, 
unlike (27). To be on safe ground, an analysis on the existence, convergence and regularity 
properties [40] of such a solution should also accompany our result. These are aspects, how- 
ever, that we will not deal with in the present work. 



Substituting (10), (19) and (20) in (21), we get 

1 



>4 



8 ( D k 

+ ^[/(M)F= T 



d 







x k 1 dx I x k 1 dx 



[/(*,*)]* 



(35) 
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whose general solution we are seeking. This non-linear partial differential equation can be 
easily reduced to a linear one, if we rewrite it as 

1 9 Wr+^M = ^l^l-lil-d/tMr+^tnl pel 



t k ~ l dt 1 k 2 I x k ~ 1 dx 1 x k_1 dx k 2 

and set 



which results in 



h(x,t) := [f(x,t)] k + ^t k (37) 



' h(x,t) = ^l\ i -^J 1 lh(x,t)} (38) 



t fc_1 dt ' <9x I x k ~ 1 



This is a parabolic differential equation. If ^- — 1 > 0, then (38) has some resemblance 
to the diffusion equation on a fractal [35], with two additional complications: the diffusion 
constant on this fractal would be time dependent 



x k-l 



(EL - l) 



t'- 1 (39) 



and there is an overall factor preceding all derivatives in the right hand side of (38). 

When D k < k, depending on the value of k, we may be able to interpret (38) as a 
time-reversed version of the diffusion equation as is frequently done in cases of dissipative 
equations. We will assume from now on, for concreteness, that D k > k. To simplify the 

subsequent equations, we set D — ^ 1 in (38). As mentioned in the first paragraph 

of this Section, to solve (38) we will follow the somewhat unusual path of using separation 
of variables instead of the more widely used Fourier transforms. We do so, in order to 
sidestep some, yet unresolved, difficulties associated to the uniqueness of the inverse Fourier 
transform [37] - [39]. Following the standard practice of separation of variables, we assume 
that 

h(x,t) = 0i(x)0 2 (t) (40) 

Substituting in (38), we find 

1 1 dxhit) 1 1 d I D dfcjx) ■ 



t k 1 02(0 dt ' 4>i(x) x k ~ 1 dx I x k ~ l dx 

where /x G 1R is a constant whose specific value is determined by the boundary conditions. 
The time-dependent part of (41) can be easily integrated and gives 

Mt) = 02(1) e^- 1 ) (42) 
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This shows that the time dependence of the solution does not grow uncontrollably as t — > oo, 
when /j, > 1 G 7Z, regardless of the value of the non-extensive parameter k. Generally 
k G R, but for most purposes of interest k > G R [2], which corresponds to k > 1 & TZ. 
The spatial part of (41) becomes 

d / 1 dM*)\ e^-Vr(x) = (43) 



[ x fc_1 dx j f) 

We set 



in terms of which (43) gives 

du{x) 
dx 



+ |Kx)] 2 -|}^- 1 = (45) 



This is a special case of the Riccati differential equation. Although it is not known how to 
solve, in general, the Riccati differential equation, (45) is, luckily, separable. 

Someone could be wondering about any deeper reasons that might justify the appearance 
of the Riccati equation (45). In geometric terms the origin of (45) in the present context 
can be traced to the following: First, we observe that (45) may be re-expressed as 

^ + u i {z) _t = o (46) 
dz D 

k i 

in terms of 1Z 3 x h-> z = G R. This is the familiar form of a Riccati equation 
that is satisfied by the logarithmic derivative of the Jacobian of a transformation [41], [42]. 
In the present case u(x) is almost the logarithmic derivative of <j>i as seen on (44). The 
presence of the multiplicative extra term x l ~ k in (44) can be traced back to the appearance 
of the term in the definition (16) which is, in turn, a result of the isomorphism (12). From 
this viewpoint, the substitution (44) becomes less mysterious. Similar things can be can 
be stated about using the g-logarithmic map (12) in going from (45) to (46). What we are 
essential doing is looking at the geometry of 1Z x 1Z from the viewpoint of R x R. The 
intrinsic geometry of each one of them is Euclidean. However the geometry of 1Z x 1Z seen 
from the viewpoint of R x R is non-trivial due to the non-linear nature of the isomorphism 
w (12). Hence the constant time lines in 1Z x 1Z do not have as inverse images constant 
time hypersurfaces in R x R. Consider the projection of a constant time line of IZxIZ 
onto a constant time line of R x R. This projection is non-trivial, as was observed above, 
and such non-triviality is expressed via the Riccati equation (45). This also explains why 
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(45) is a Riccati equation with respect to the tangent vector to the constant time line, rather 
than being a differential equation with respect to time (normal to the constant time line) 
direction itself as in [41], [42]. 

We solve first the special case of (45), for which 

u\x) = t (47) 
for which (45) actually becomes linear. Then 

i 



f(x,t) = e^*- 1 ) - — t k j (48) 

Here A(fi) and x Q are constants to be determined from the initial conditions. Due to the 
linearity of (45), and the lack of any additional constraints for /i, the general solution of 
(45) is formally given by 

oo 

F(x,t) = J f(x,t) <£> <V (49) 
i 

In order for F(x, t) to be interpreted as a probability distribution function it should satisfy 

F(x, t)>l, V (x, t) E U x 11+ (50) 
and it should be normalized in TZ, namely it should satisfy 

+oo 

J F(x,t) ®d®x = + (51) 

— oo 

In case F(x, t) is not automatically normalized, the normalization can be achieved by 
setting 



p(x,t) = F(x,t) N (52) 



where 



+oo 

N = 



In the general case when 
then (45) becomes 



J F(x,t) I) d®x (53) 



u 2 (x) (54) 



D 



du A ._j 



, = dx (55) 

u*(x)-% 
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which upon integration gives 



x k -x k i — nr x k -x k 



[AM + M eV%- X ~^ + [AM - M e~V% : 



ul(l x\ — * / ~~ * — — f56l 

V D [AM + e VV-*- - [AM - JI] eSs—* 

Using (44) and integrating once more, we find 



Mx) = 2 Bin) Jt — —-j-r (57) 

V D [AM + y/%]- ^M) - yj%\ e- 2 ^ ■ 

where BM is another integration constant that is determined by the initial conditions. 
Substituting (57) and (42) into (37), we eventually find 

ji . . evo k ek v >-) ft _ i k 



x k -x k 



[AM + M - [AM - M e- 2 ^ 



f(x,t) = 2^ CM — — y^-, ~ -jr- ** 

(58) 

where C(/x) = i3(/i)0 2 (l)- The general solution to the diffusion equation (21) is given by 
substituting (58) into (49) for a given set of initial conditions. The required normalization 
is, once more, provided by (52) and (53). 

We observe that (58) has the form of a stretched exponential reminiscent of solutions of 
the porous medium equation to which it may worth comparing. A second, related, reason is 
that its known solutions exhibit scaling behaviour [43], [44] which is an important property 
shared by structures described by the Tsallis entropy, such as fractals, for instance. A third 
reason for singling out and paying particular attention to solutions of the porous medium 
equation, among all the possible non-linear and fractional generalisations of the diffusion 
equation, is because the porous medium equation arises as the gradient flow of the Tsallis 
entropy functional [45]. 

Foregoing all specifics, which can be found in standard sources [46], the porous medium 
equation is the non-linear second order quasi-parabolic equation given by 

^|^ = V 2 [^(f,t)r, m>\ (59) 

The Zel'dovich-Kompaneets/Barenblatt (ZKB) solution is a radially symmetric, self-similar 
solution having the form, in M. N x IR + 

^,t)~f Ct2C " A||g||2 V" 1 (60) 
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where 

C=(Nm-N + 2)-\ X= (m ~ 1)C (61) 

2m 

and c > is a constant depending on the total mass 

M — f g(x,t) d N x (62) 

The significance of the ZKB solution (60) is considerable if one takes into account [47] that 
under initial data obeying certain conditions, the solutions of the porous medium equation 
converge to (60) for t — > oo. We clearly see two common features between (48), (58) and 
(60). The first one is the appearance of the same exponent, given in two different notations, 
in the overall solution. The second is the appearance of a time- dependent "drift" term in 
both cases. Naturally, there are differences between (48), (58) on the one hand and (60) 
on the other: one of them is that (48), (58) have the form of s stretched exponential in the 
spatial distance, whereas (60) has the form of a stretched Gaussian with respect to analogous 
variable. Moreover, the temporal asymptotic behaviour of (48), (58) is exponential, whereas 
that of (60) is power-law/polynomial. In addition, we observe that it is the logarithm of 
(48), (58) that exhibits, asymptotically, scaling behaviour as x,t — > oo. Equivaletly, it is 
the exponent of (48), (58) that scales nicely as x, t — >■ oo. Indeed, under the transformation 

x i — y ax, t i — y at (63) 

we observe that 

k 

log /CM) H- °- 2 log/OM) (64) 

This scaling behaviour can be traced back, once more, to the exponential form of the iso- 
morphism (12) relating R and 71. For k > 1 (48) and (58) can also be interpreted as 
log-convex measures on WL N , which indicates that they have common features with the 
Lebesgue and the Gaussian measures in M. N such as their concentration properties [48]. 
Measure concentration is particularly important [49], since it can be seen as one important 
reason why the predictions of Statistical Mechanics result in effectively constant values of the 
experimentally observed quantities, at least in the thermodynamic limit. The concentration 
of a measure which is not a tensor (power) of another one, allows us to extend results that 
stem from the conventional assumption of independence to far weaker inter-dependencies 
[50]. In this context, it may be worth exploring the concentration properties of measures 
induced by structures reflecting the composition properties of the Tsallis entropy (3). It 
may also be with exploring the relation of such measures, if any, to a recent generalization of 
the Central Limit Theorem [16] - [19] motivated by and tailored to the needs of the Tsallis 
entropy induced viewpoint toward independence and additivity. 
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4. Discussion and conclusions 



We have found formulated a generalized diffusion equation on the set 7Z whose def- 
inition, algebraic and differential operations were inspired by the non-extensive properties 
of the Tsallis entropy [28]. We determined its general solution which happened to be a 
stretched exponential, asymptotically, in both the spatial and time variables. We compared 
the asymptotic form of our solution to that of the ZKB solution of the porous medium 
equation. We did all this in order to get a better understanding of the implications of the 
Tsallis entropy composition property (3) and subsequently of (4) as compared to the the 
BGS composition of independent events, which is the regular addition. 

An extension of this work is motivated by the realisation that there is not really any 
compelling reason why the temporal and spatial behavior of a system should be described 
by the same value of the deformation parameter k in (21), at the mesoscopic level. If the 
evolution of a system is described by a Riemannian or Finslerian metric on its phase space, 
and if such a metric has a symmetry group that reflects the treatment of the spatial and tem- 
poral directions on the same footing and with the same parameters, it may be advantageous 
to preserve these symmetry properties in the system's coarse grained/statistical description. 
It is important to recall that the diffusion equations are inherently non-relativistic, so they 
manifestly break the, potential underlying Lorentz invariance, if such invariance even exists, 
of the microscopic system. For this reason, it may not be unreasonable to allow for different 
values of k in the spatial and temporal directions of mesoscopic equations. Such equa- 
tions may arise more naturally from the composition properties of entropic functionals that 
depend on two parameters [2]. In such cases, an obvious generalization of (21) would be 

k 

d m f(x,t) = D <g) d <S)X d <S)X f(x, t) (65) 

with K 7^ k. It might be of some interest to determine the general solution of (65), if 
this is practically feasible at all, and ask for its asymptotic behavior as well as its physical 
significance, if any. Otherwise, we may want to seek solutions of (65) by using the "Tsallis 
entropy- inspired" q-exponential ansatze [2]. 

As a result of the above analysis, we start seeing a few points that may be more generic 
and of broader applicability that supersedes the way they were reached in this work. When 
compared to the ordinary diffusion in M. N , the solutions (48), (58) grow much faster in time. 
Hence they are more akin to super- diffusion rather than to regular diffusion. If we still 
wanted to emulate such a process by a diffusion, we would say that it would have to take 
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place not on R but on a space that was negatively curved. This since on such a space 
any quantity would diffuse faster than in M. N , as should be intuitively obvious [35]. Hence, 
form the viewpoint of ordinary diffusion, it is as if the use of 1Z amounts to endowing M. N 
with a negative curvature metric. So, we arrive at a previously reached result [31], where the 
Tsallis entropy composition (3) induces a hyperbolic structure on the underlying (topolog- 
ical) space. It seems therefore that the underlying hyperbolicity encoding (3) is somewhat 
insensitive to the choice of the base field, be it 1Z as in the present case or M. q as in [31]. 

A far more general application of this diffusion approach with important implications for 
geometry be can be followed. One could reverse the above line of reasoning: since the diffu- 
sion properties ultimately depend only on the Ricci and not on the sectional curvature per 
se, one could use the behaviour of the Tsallis entropy, and especially its convexity [45], to de- 
fine the Ricci curvature for non-Riemannian spaces [51]. One such class of non-Riemannian 
spaces that has attracted considerable interest relatively recently in Statistical Mechanics, 
are networks/graphs [52] - [55]. A recently defined Ricci curvature [56] for discrete metric 
spaces, such as graphs, provides an analytical tool for potentially unraveling their pertinent 
properties. This can be particularly useful especially in the research direction of networks 
in which only a very small subset of its results have been analytically obtained [54], [55] and 
illustrates the robustness as well as the wide range of applicability of the Tsallis entropy 
formalism. 
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